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Abstract — In this letter, we study the robust beamforming 
problem for the multi-antenna wireless broadcasting system with 
simultaneous information and power transmission, under the 
assumption of imperfect channel state information (CSI) at the 
transmitter. Following the worst-case deterministic model, our 
objective is to maximize the worst-case harvested energy for 
the energy receiver while guaranteeing that the rate for the 
information receiver is above a threshold for all possible channel 
realizations. Such problem is nonconvex with infinite number of 
constraints. Using certain transformation techniques, we convert 
this problem into a relaxed semideflnite programming problem 
(SDP) which can be solved efficiently. We further show that the 
solution of the relaxed SDP problem is always rank-one. This 
indicates that the relaxation is tight and we can get the optimal 
solution for the original problem. Simulation results are presented 
to validate the effectiveness of the proposed algorithm. 

Index Terms — Energy harvesting, beamforming, worst-case 
robust design, semidefinite programming. 



I. Introduction 

Energy harvesting for wireless communication is able to 
extend the flying power of handheld devices and advocacy for 
green communication |[T|- Q. With the aid of this promising 
technique, the transmitter can transfer power to terminals 
who need to harvest energy to charge their devices, which 
is especially important for energy-constrained wireless net- 
works. Beamforming is another promising technique which 
exploits channel state information (CSI) at the transmitter 
for information transmission ID- Q. In wireless networks 
with simultaneous transmission of power and information, 
beamforming is anticipated to play an important role as well. 

The beamforming design with perfect knowledge of CSI 
at the transmitter was first considered in Q to characterize 
the rate-energy region in a simplified three-node wireless 
broadcasting system. In practical scenarios, perfect knowledge 
of CSI may not be available due to many factors such as 
inaccurate channel estimation, quantization error, and time 
delay of the feedback. 

The goal of this letter is to investigate the robust beam- 
former design with imperfect CSI for simultaneous informa- 
tion transmission and energy harvesting. In general, there 
are two classes of models to characterize imperfect CSI: 
the stochastic and deterministic (or worst-case) models. In 
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the stochastic model, the CSI errors are often modeled as 
Gaussian random variables and the system design is then based 
on optimizing the average or outage performance H), ||9|. 
Alternatively, the deterministic model assumes that the CSI 
uncertainty, though not exactly known, is bounded by possible 
values ifTOl . lITTl . In this case, the system is optimized to 
achieve a given quality of service (QoS) for every possible 
CSI error if the problem is feasible, thereby, achieving absolute 
robustness. It was also shown in lfT2l that a bounded worst- 
case model is able to cope with quantization errors in CSI. In 
this letter, we shall employ the worst-case approach to address 
the robust beamforming design problem. 

Consider the three-node system shown in Fig. 1, where 
we assume that the transmitter only has imperfect knowledge 
of the channels to both the information receiver and energy 
receiver We formulate the worst-case robust beamforming 
problem for harvested energy maximization at the energy re- 
ceiver while ensuring a minimum target rate at the information 
receiver Since the original problem has infinite constraints due 
to the channel uncertainties, we first transform it into an easier 
problem which has finite constraints but is still nonconvex. 
Then we apply the semidefinite relaxation (SDR) and obtain 
a semidefinite programming (SDP) problem which can be 
solved efficiently. Finally we show that the optimal solution 
of the SDP problem is always rank-one, which means that the 
relaxation is tight and we can obtain the optimal solution of 
the original problem. 

The rest of this letter is organized as follows. In Section II, 
the system model and the problem formulation are presented. 
Section III presents our proposed algorithm to find the so- 
lutions to the robust problems using convex optimization and 
rank relaxation, and show its optimality. Simulation results are 
given in Section IV. Finally, Section V concludes this letter. 

Notation: (•)^ and Tr{ } stand for Hermitian transpose and 
the trace respectively. \x\ denotes the absolute value of the 
scalar x and ||x|| denotes the EucUdean norm of the vector x. 
The function log(.) is taken to the base 2. 

II. System Model and Problem Formulation 

With reference to Fig. 1, we consider a three-node multiple- 
input single-output (MISO) communication system, where the 
transmitter has N antennas and each receiver has a single 
antenna. Let and denote the frequency-flat quasi-static 
1 X complex channel vectors from the transmitter to the 
information receiver and the energy receiver respectively, and 
s denote the transmitted symbol. Then the received signals at 
the information receiver and the energy receiver are given by. 
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Energy Receiver 



Information Receiver 

Fig. 1. A three-node wireless MISO system for simultaneous information 
and power transmission. 
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g^ws + Zg 



(1) 
(2) 



where w is the x 1 beamforming vector appUed to the 
transmitter, and and Zg are the additive white circularly 
symmetric Gaussian complex noise with variance on 
each of their real and imaginary components. 

For the energy receiver, it will harvest energy from its 
received signal. Thanks to the law of energy conservation, we 
can assume that the total harvested RF-band power, denoted 
by Q, is proportional to the power of the received baseband 
signal, i.e., 

= 77|g«w|2 (3) 

where 77 is the efficiency ratio at the energy receiver for 
converting the harvested energy to electrical energy to be 
stored. Here we simply assume that rj = 1 and the details 
for the converting process is beyond the scope of this letter. 

Our objective is to maximize the harvested energy for the 
energy receiver while guaranteeing that the information rate 
for the information receiver is above a threshold. Mathemati- 
cally, the problem is expressed as follows: 



\ H |2 

max |g w| 



s.t. log 1 + 



> r 



< P 



(4) 

(5) 

(6) 



where r is the rate target for the information receiver and 
P is the power constraint at the transmitter. Similar problem 
has been considered in |i7J with the objective of maximizing 
information rate subject to a minimum energy threshold. 

Herein we consider that the transmitter has imperfect CSI 
of both receivers. In particular, the channels are modeled as 



h = h + Ah 

g = g + Ag 



(7) 
(8) 



where h and g denote the estimated CSI known at the 
transmitter. Ah and Ag are the error vectors. We assume 
no statistical knowledge about the error vectors but that they 
are bounded by some possible values (also known to the 
transmitter) as 



|Ah||<£ 

|Ag|| < e 



(9) 
(10) 



where e is the radius of the uncertainty region. We assume 
that both the two receivers have perfect CSI knowledge. 

To take the CSI errors into account, the problem Pq based 
on worst-case criterion can be formulated as 



Pi : max min |(g + Ag)^wP 

w ||Ag||<e 



(11) 



s.t. log(l + K£±^^)>.,Vi|Ah||<s(12) 



< P. 



(13) 



Since log(l + x) is monotonically increasing for positive x, 
problem Pi can be reformulated as below 
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Pi : max min |(g + Ag)^w| 

w ||Ag||<e 



(14) 



s.t. |(h + Ah)^w|2 > (72(2'^ - l),V|lAh|| < e (15) 



< P. 



(16) 



It can be seen that the goal of the problem Pi is to maximize 
the harvested energy for the worst channel realization while 
guaranteeing that the information rate is above a threshold for 
all possible channel realizations. 

III. Semidefinite Programming Solution 

The key challenges in problem Pi are the channel un- 
certainties and the nonconvex constraints, which cause that 
Pi is a semi-infinite nonconvex quadratically constrained 
quadratic programming (QCQP) problem. It is well known 
that the general nonconvex QCQP problem is NP-hard and 
thus, intractable. However, as we will show in the following, 
due to the special structure of the objective function and the 
constraints, problem Pi can be reformulated as a convex SDP 
problem and solved optimally. 

We first transform the above problem into a more tractable 
form. For the objective function of Pi in (fT4l i. we simplify 
it using an approach similar to the one developed in flOl and 
ifTSl . According to triangle inequality, we obtain 



|g^w + Ag^wl > |g^w| - |Ag^w|. 



(17) 



Then applying the Cauchy-Schwarz inequality to the second 
term in the right-hand-side (RHS) of (fTTI i. we have 

|Ag^w|<!|Agi|.||w||<ej|wj|. (18) 

Plugging ( fTSl l into (fTTI i. we then have that 

|g^w+Ag^w| > |g^w|-|Ag^w| > |g^w|-e|lw|l. (19) 

An important observation about problem Pi is that its optimal 
solution is obtained only when the constraint in (fT6] l is active, 
i.e., the transmitter should work with full power Then we have 



Ag« 



w| > 



li'^wi 



:\/p. 



(20) 



The inequality becomes equality when Ag = ^ ^^^e 
where 6 is the angle between g^ and w. Note that it has been 
assumed that |g'^w| > £||w|| in ( fT9] l, and |g^w| > e^/P in 
( [20I 1. This assumption essentially means that the errors Ag is 
sufficiently small or equivalently e is sufficiently small. It is 
a practical assumption since large channel estimation errors 
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can cause large beamforming errors and no robustness can 
be guaranteed in such case. Then combining (fT7]i-([20]i. we 
conclude that 

2 

min |(g + Ag)^w|2= \g"w\-eVP . (21) 

l|Ag||<e 

For the infinite number of constraints in ( fTSl l. we can 
similarly have that 

|h^w + Ah^w| > Ih-^wl - eVP. (22) 

Here, the equality holds when Ah = with ip 

being the angle between and w. Then in order to meet 
the constraints for all possible Ah, we just need to satisfy the 
following 

|h^w| - sVp > o^j2^ - 1. (23) 

Then the robust beamforming problem Pi can be rewritten 
as follows 

Pi : max |g^wp 



t. Ih^wp > ( e/P + (tV2'' - 1 



llwf < P- 

Although the problem Pi is much easier now, it is still a 
nonconvex QCQP problem. We then apply the semidefinite 
relaxation and obtain the following relaxed problem; 



P2 : max Tr{GW} 

w>-o 



(27) 



s.t. Tr{HW} > (^eVP + (tV2'- - ij (28) 

Tr{W} < P (29) 

where G = gg^ and H = hh^. Notice that the rank-one 
constraint has been dropped and P2 is a relaxed version of Pi. 
The problem P2 is a standard SDP problem which is convex 
and can be solved efficiently using the software package lfT4l . 

At this point, an important question is that whether the 
optimal solution of P2 is rank-one. If W is rank-one, then 
the optimal beamformer for the original problem Pi can 
be extracted by eigenvalue decomposition. Otherwise, the 
solution of P2 is only an upper bound of Pi and the 
beamformer extracted from W is not guaranteed to be globally 
optimal. Generally there is no guarantee that an algorithm 
for solving SDP problems will give the desired rank-one 
solution. However, in some special cases such as ifTSl - ifTTl , 
the relaxation is proven to be exact and thus there always 
exists a rank-one solution. Whether the relaxation is tight for 
our proposed algorithm will be addressed in the following 
theorem. 

Theorem 1: The optimal solution W for problem P2 is 
rank-one. 

Proof: Please refer to Appendix A. ■ 
According to Theorem 1, we can see that problem P2 is 
indeed equivalent to the original problem Pi, which means 
that the relaxation is tight. So in order to solve the problem 
Pi, we first solve the SDP problem P2 and obtain the resulting 
rank-one matrix W*. Apply eigenvalue decomposition on W* 
as 
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Fig. 2. Average harvested energy for the robust beamforming design. 
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Fig. 3. Outage percentage for the nonrobust beamforming design. 



The optimal solution of Pi is then obtained as w = \/a*w*. 

IV. Simulation Results 

In this section, we present numerical results to evaluate 
the performance of the proposed robust beamforming algo- 
rithm. We consider the three-node MISO system in which the 
transmitter has four antennas (N = 4). We set the power 
P = 10 and noise covariance <t^ = 1. The channel from 
the transmitter to each receiver is assumed as the normalized 
Rayleigh fading channel. The rate target for the information 
receiver is set to be smaller than log(l + 10||h|p) in order 
for the problem to be feasible. For simplicity, we normalize 
these channel vectors with respect to the number of transmit 
antennas as ||h|p = ||g|p = 4. Thus the feasible region 
for the rate target is < r < log(l + 40). A total of 100 
independent normalized channel realizations are simulated. 
For each channel reaUzation, 100 channel uncertainty samples 
are generated. 

In Fig. 2, we plot the average harvested energy versus differ- 
ent targets of information rate for different levels of bounded 
channel uncertainty. The special case with the CSI (i.e., e = 0) 
is also simulated. It can be seen that the performance loss is 
small when the CSI error is not big. Also the performance gap 
increases when the rate target becomes larger 

In order to show how important it is to take the channel 
uncertainty into account when designing the beamformers, let 
us assume that the beamforming design takes place under the 
assumption of perfect CSI at the transmitter while in fact 
there is some uncertainty associated with the CSI used in the 
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design problem, which we call the "nonrobust beamforming 
design". Fig. 3 shows percentage of outage Q at different rate 
targets for the nonrobust design. We observe that the channel 
uncertainty, when not considered in the design process, leads 
to frequent violations of the rate target at the information 
receiver However, for our proposed worst-case robust beam- 
forming algorithm, the rate target is always satisfied and no 
outage happens. 

V. Conclusion 

In this letter, we consider the worst-case robust beamform- 
ing design for the wireless communication system with both 
information and energy receivers when the CSI is imperfect. 
By means of semidefinite relaxation, we transform the original 
robust design problem into a SDP problem. Then we prove 
that such relaxation is tight and we can always obtain the 
optimal solution of the original problem. The performance of 
the proposed beamforming algorithm has been demonstrated 
by simulations. Future research directions may include the 
robust beamforming design for the more general broadcast- 
ing systems with multiple information receivers and multiple 
energy receivers. 

Appendix A 
Proof of Theorem 1 

Denote l3 = (^e\^ + ay/T^~fj , the Lagrangian of P2 
is given by 



L (W, A, At) = Tr{GW} + A (^Tr{HW} - /3J (31) 

-/.(Tr{W}-P), 

where A and fi are the dual variables. The Lagrange dual 
function is then defined as 



max L (W, A, 11) . 

W)-0 



(32) 



Since P2 is convex with strong duality, we can solve it by 
solving its dual problem 



D2 : mill q(X, u). 

A>0,Ai>0' 



(33) 



Denote the optimal solution of D2 as {X*,fi*), then the 
matrix W* that maximizes L(W, A*,/i*) is the optimal 
solution of P2, which means that we can find W* through 
the following problem 

max Tr{GW} - Tr | - A*h) w| (34) 

where the constant terms has been discarded. In order for 
problem (|34] | to have a bounded value, it is shown as follows 
that the matrix — A*H should be positive definite. Suppose 
1^1*1 — A*H is not positive definite, then we can choose W = 
tww", where i > and TrJt{p.*I- X*'H.)ww^} < 0. Due to 
the independence of g and h, it follows that Tr{<Gww^} > 
with probability one. Let t +00, the optimal value in ( [34l i 
will be unbounded, which is a contradiction of the optimality 

of ix\^l*). 

'We call the outage happens when the rate target is not satisfied at the 
information receiver. 



Define Q = {fi*I - A^H) and let W = Qi/^WQi/^ 
the problem in (l34l i is then rewritten as 

-i/ss^iffwi-n-i/ag^ _ Tr{W}. 



max (Q"'/'g)''W(Q- 
w>-o 



(35) 



Then we claim that the optimal solution of ( |35] ) is always 
rank-one. Suppose the optimal solution W is not rank- 
one, without loss of generality, we can assume its rank is k 
{2 < k < N) and decompose it as = '^f^j^f ■ 

Then we choose another W = (X]j=i c>^i)'^i'^f^' where 

i = arg max |(Q~^^^g)'^Wj |. Then W can achieve a 

ie{i,...,fc} 

larger value than W , which is a contradiction. 

From the above discussions, it is known that W is always 
rank-one. Since W* = Q~^/^W Q~^/^, we must have that 
W* is rank-one, which completes the proof of Theorem 1. 
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